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Abstract 

The time equation associated to the Dirac Equation (DE) is stud- 
ied for the radiation-dominated Friedmann-Robertson- Walker (FRW) 
universe. The results are analysed for small and large values of time. 
We also incorporate the corrections of the paper studied by Zecca [1] 
for the matter-dominated FRW universe. 
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1 Introduction 



It is of cosmological interest to consider a universe of incoherent radiation. 
The primordial radiation played a dominant role in the early universe for the 
times smaller than the times of re-combination. Zecca [1] has discussed the 
time equation associated to the DE in the FRW spacetime in the case of Stan- 
dard Cosmology. He confined his paper to the study of matter-dominated 
cosmological model only. This paper is dedicated to a further study of the 
time equation associated to the DE of a radiation-dominated cosmological 
model. Also, we will incorporate the corrections in the expansion terms of 
Eqs.(14) and (17) in the Zecca's paper [1]. 

It is well known that the DE can be written in General Relativity in terms 
of covariant derivatives and generalised Pauli matrices. In the context of the 
Newmann- Penrose formalism [2], the DE can further be expressed in terms 
of directional derivatives and spin coefficients. The DE, formulated by the 
spinorial formalism of Newmann and Penrose [3] , is given by 



particle [4]. If we correspond <-> (P, Q), ip <-> (U,V), we can define the 
spinor 



define inner product between the solutions of the DE by setting [5] 
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The Eq.(3) is independent of the spacelike Cauchy hypersurface £ of volume 
element oE, n a is a future directed unit vector orthogonal to dT,, <j aa , is 

the generalised Pauli matrix and has the form a AA , = \ 

The separation of the angular part in Eq.(l) is performed by a standard 
separation method which gives the Teukolski-like equation for spin(l/2) field 
admitting explicit solution [6]. The surviving coupled equations in the time 
and radial variables are separated by writing the unknown wave function in 
terms of a known particular solution. These equations can be expressed in 
the conformal time parameter r and in the spatial parameter s defined by 

The complete spinorial solution of Eq.(l) has been found and is of the form 

pi = $^H(r,t){(-iyA lk (s)T k (r) + [fA lk (s) - A' lk (s)} drT k }, ) 
Qi = ~I^H(r,t){A k (s)T k (r) + (-iy[fA lk (s) - A' lk (s)} £ dfT k }, J 

(7) 

where i — 0,1. The angular functions appearing in Eq.(7) are of the form 
Si£l = S\H(9 , 4>)' 's(j = 1,2), where Si^(0,<f>) are solutions of an eigen-value 
problem originated by the solution of the angular part of Eq.(l). The Ai k (s)'s 
are the solutions of the angular and radial equations, respectively whose 
explicit expressions can be found in [7]; \,k are separation constants such 
that A 2 = (/ + l) 2 , 1 = 0,1, 2, ... for m = 0; A 2 = (I + 1/2) 2 , 1 = \m\, \m\ + 
1, \m\ + 2, ... for m — 0, ±1, ±2, ±3, .... The angular functions are assumed 
to satisfy the normalization condition 

J dnsiiiio, 4>)s\P m ,{e, <j>) = 5 lv 5 mm , (8) 

The function H is connected to a particular integral of Eq.(l) and has the 
form 




H(r,t) 



1 fl + VI - ar 2 



R l ' 2 {t) 



exp[in*V2t], (a = 0,±l). (9) 



The function Tfc(r) is a solution of the separated time equation in the 
conformal time parameter r and satisfies [7] the equation 

T" + 2V2l^R{t)T' + {2lV2^R'{t) + k 2 )T = 0, (10) 

where prime denotes differentiation with respect to r. The solution of Eq.(10) 
depends on the dynamical evolution of the cosmological background and this 
satisfies a property similar to the Wronskian property of the scalar field time 
equation in FRW spacetime [8,9]. By using a first formal integration of 
Eq.(10), one can easily show that any solution Tkij) of Eq.(10) satisfies the 
relation 

\T k (r)\ 2 + k 2 

This normalization is a necessary requirement of a second quantization of the 
Dirac field. It is noted that Eq.(10) can be solved in particular simple cases 
such as that of static universe or of the neutrino case. However, it seems 
difficult to give the solution of Eq.(10) for a general cosmological evolution. 

The layout of this paper is as follows. In section two, we analyse the 
matter-dominated FRW universe to incorporate the corrections studied by 
Zecca [1]. In the next section we shall extend this procedure to study the 
solutions for the radiation-filled cosmological model. Section four is devoted 
to discussion of the results. 

2 Solutions in the Matter-Dominated Cos- 
mological Model 

It is well known that the FRW spacetime whose metric 



T k (?)dt 



= constant = 1. 



(11) 



ds 2 = dt 2 - R 2 {t) 



dr 1 



1 — ar' 



+ r 2 (^ 2 + sin 2 0d<p 2 ) 



(a = 0,±l) (12) 



is the base of the standard model [10,11]. This is the natural context for 
the formulation of the Einstein Field Equation for a general isotropic and 
homogeneous universe. We shall first consider the matter-dominated universe 
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for the purpose of correction in [1]. For matter-dominated universe, R(t) in 
parametric form is given as 



R(t) = Et 2 '\ (a = 0), 

R(t) = A(cosh^- 1), t = B(smhtp-ip), ip > 0, (a = -1), I (13) 
R(t) =C(l- cos 9), t = D(e-am6), 0<6<2n, (a = 1). 

In terms of the conformal time parameter r defined in Eq.(6) we obtain 

R(r) = £ 3 r 2 /9, (a = 0), 1 

i?(r) = A[cosh(TA/ J B) - 1], (a = -1), I (14) 

R(t) = C[l - cos(rC/D)], (a = 1). J 

If we use Eq.(14) in Eq.(10) they all give a linear differential equation with 
finite regular point of the form 

T" + ( £ p n A T' + ( £) g„r") T = 0. (15) 

\n=0 / \n=0 / 

By using a standard method [12], we can set 

oo 

T = £ cnr", (16) 

?1=0 

where the coefficients c n can be found by the recurrence relation 



(n + 2)(n + l)c n+2 + £^(n + 1 - j) c „+i-j + £ ^ c «-i = °> n = 0, 1, 2, ... 

j'=o i=o 

(17) 

We shall always take c\ — due to our interest only in the time integral. 
For the flat model a = 0, the recurrence relation (14) in [1] will become 



T = Co {l - ^k 2 - ^2aE 3 /9 + ^k 4 + ^Uak 2 E 3 /9 

+£{40(aE 3 /9) 2 - A; 6 ] - ^44aE 3 A; 4 /9 + ^[A; 8 - 628A; 2 (a£ 3 /9) 2 ] 
+ ^[100a£ 3 A; 6 /9 - 2240(a£ 3 /9) 3 ] + ...}, 

(18) 

where a = 2\/2 i/a*. 
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For the open model a = — 1, the recurrence relation (17) in [1] will take 
the form 



T = Co {l - ^k 2 - faA 3 /B 2 + ^k A + ^[7aA 3 k 2 /B 2 - aA 5 /B 4 ] 
+ £{lO(aA 3 /B 2 ) 2 - k 6 ] + ^(16aA 5 k 2 /B 4 - 22uA 3 k^B 2 - aA 7 / B 6 ) 
+^[-157(aA 3 /B 2 ) 2 k 2 + 5Qa 2 A s /B 6 + A; 8 ] + ...}. 

(19) 

By combining Eqs.(10),(13) and (14) one has T = R- 1 = t' 2 / 3 = r~ 2 for 
t>1. Thus by using Eqs.(7),(9) and (18), the time dependence of H will 
take the form 

r 1/3 exp[L^V2t], (t > 1, a = 0). (20) 



3 Solutions in Radiation-Dominated Cosmo- 
logical Model 

In this section we shall extend the above procedure to solve Eq.(10) for 
radiation-filled FRW universe model. For radiation-dominated FRW space- 
times, the scale factor R is given as 



R(t) = Et 1 / 2 , (a = 0), 

R(t) = Asinh^, t = B(coship - 1), tp > 0, (a = -1), 
R(t) = Csm0, t = D(l -cos0), < 9 < 2n, (a = 1). 



(21) 



We can write these values of scale factor in terms of the conformal time 
parameter r as 



R(t) = E 2 t/2, (a = 0), 
R(t) = Asmh(rA/B), (a = -1), 
R(t) = Csm(rC/D), (a = 1). 

(i) The Flat Model a=0 

By using Eqs.(10), (15) and the 1st equation in (22), we can have 

Po = 0, pi = aE 2 /2, p n = 0, n ^ 1, 1 
q = aE 2 /2 + k 2 , q n = 0, n > 1. J 



(22) 



(23) 
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With the help of Eq.(23) in the recurrence relation (17), we can find all values 
of c. Replacing these values of the c's in Eq.(16), we get 

T = Co {l - ^(aE 2 /2 + k 2 ) + T -^(3aE 2 /2 + k 2 )(aE 2 /2 + k 2 ) 
-£{5aE 2 /2 + k 2 ){3aE 2 /2 + k 2 ){aE 2 /2 + k 2 ) 
+ ^(7aE 2 /2 + k 2 )(5aE 2 /2 + k 2 )(3aE 2 /2 + k 2 )(aE 2 /2 + k 2 ) + ...}. 

(24) 

(ii) The Open Model a=-l 

From Eqs.(lO), (15) and the second in Eq.(22) with the use of the series 
expansions of cosh a;, sinhx, we can make the identifications 



Po = 0, p 2n . 1 = j^ w (A/B) 2n -\ n>l, 

q = aA 2 /B + k 2 , q 2n = ^(A/B) 2 \ n > 1. 




Making use of these values in the recurrence relation (17) to find c's and then 
Eq.(16) yields 

T = Co {l - ^(aA 2 /B + k 2 ) + ^[(3aA 2 /B + k 2 )(aA 2 /B + k 2 ) 

-a{A 2 /B){A/B) 2 ] - ^[{5aA 2 /B + k 2 ){{3aA 2 /B + k 2 ) 
(aA 2 /B + k 2 ) - a{A 2 /B){A/B) 2 } - l0a(A 2 / B)(A/ B) 2 
(aA 2 /B + k 2 ) + a(A 2 /B)(A/B) 4 ] + ^[(7aA 2 /B + k 2 ) 
{(5aA 2 /B + k 2 ){(3aA 2 /B + k 2 )(aA 2 '/B + k 2 ) (26) 
-a(A 2 /B)(A/B) 2 )} - 10a(A 2 /B)(A/B) 2 (aA 2 /B + k 2 ) 
+a(A 2 /B){A/Bf} - 35a(A 2 /B)(A/B) 2 {(3aA 2 /B + k 2 ) 
(aA 2 /B + k 2 ) - a{A 2 /B){A/B) 2 )} + 3a{A 2 / B){A/ Bf 
(aA 2 /B + k 2 ) - a{A 2 /B){A/Bf] + ...}. 

(iii) The Closed Model a=l 

The closed model can be obtained from the open model (a = —1) by 
replacing A by iC and B by —D. Using the normalization Eq.(ll), we can 
have c = 1 in every case. 
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4 Discussion of the Reults 



We see from Eqs.(24) and (26) that when r tends to zero T(r) approaches 
to 1. This implies that the complete spinorial solution (7) becomes inde- 
pendent of the conformal time parameter r for r — > 0. When we combine 
Eqs.(7),(9),(21),(22),(24) and (26), we obtain the same behaviour in every 
case for t « 1. 

Consider the case of a cosmological model with a = 0. For large t,R = t 1 / 2 
hence R = tso that the asymptotic behaviour of the solutions of Eq.(10) are 
now approximated by T = r^ 1 = t~ 1 / 2 so that H has the time dependence 

r 1/A exp[ifi*V2t}, (t > 1, a = 0) (27) 

Similarly, for the open (a = —1) radiation-dominated universe we have R = t 
and hence R = exp(rA/B) for large t. Thus the asymptotic behaviour of the 
solutions of Eq.(10) is given by T = exp(—rA/B) = t^ 1 and for H the time 
dependence is of the form 

r 1/2 exp[in*V2i\, (t > 1, a = -1) (28) 
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